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Abstract. We construct a family of triatomic models for heteronuclear and 
homonuclear molecular Bose-Einstein condensates. We show that these new 
generalized models are exactly solvable through the algebraic Bethe ansatz method 
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1. Introduction 

Since the pivotal experimental achievement that led to realizations of Bose-Einstein 
condensates (BECs), using ultracold dilute alkali gases [H [2], a great effort has been 
devoted to the understanding of new properties of BEC. The creation of a molecular 
BEC compound, which has been obtained by different techniques [3], leads one also to 
the compelling chemistry of BECs, where the atomic constituents may form molecules, 
for instance, by Feshbach resonances [I] or photoassociation [5|. These results turned 
the search for integrable models that could be candidates for describing BEC properties 
into a very active field of research [6J [TJ, El El EDI EH EG]. In fact, exactly solvable models 
are expected to provide a significant impact in this area, a view that has been promoted 
in [T31CE1]. 

Moreover, the recent experimental evidence for Efimov states in an ultracold cesium 
gas [15] provided a physical ground for the search of tri-atomic homonuclear molecular 
BECs. Due to the rapid technological developments in the field of ultracold systems, 
this experiment is just the beginning of the study of tri-atomic molecules [HI ITT] . 

We introduce, in the present paper, a complete family of new solvable models 
for both heteronuclear and homonuclear tri-atomic molecular BECs obtained through a 
combination of two Lax operators constructed using special realizations of the su(2) and 
su(l, 1) algebras, as well as a multibosonic representation of si (2), discussed recently in 

Until now, integrability was shown for those Hamiltonians describing heteronuclear 
and homonuclear di-atomic molecular BECs [Tj. In the heteronuclear case, two different 
atoms, labelled a and b, can be combined to produce a molecule labelled by c. The 
different degrees of freedom in such models are represented by canonical creation 
and annihilation operators {a, b, c,...,a*, b\ c' , ...} satisfying the usual commutation 
relations [a, at] = I, etc., and commuting among different species. The Hamiltonian for 
the di-atomic heteronuclear model reads [TJ [19] 

H = U aa N 2 a + U bb N 2 b + U CC N 2 C + U ab N a N b + U ac N a N c + U bc N b N c 

+ fi a N a + fi b N b + /j c N c + n(a)b ] c + c ] ba). (1) 

The parameters describe S-wave scattering, ^ are external potentials and fl is the 
amplitude for interconversion of atoms and molecules. In the homonuclear case, two 
identical atoms labelled a are combined to produce a molecule c [TJ . The generalization 
to the triatomic case is not immediate and in particular, for the homonuclear case, a 
one mode multibosonic realization of sl{2) is essential. 

In what follows we will show that it is possible to find other integrable generalized 
molecular BEC models, more specifically, heteronuclear and homonuclear tri-atomic 
models. We first introduce two models, one with two identical species of atoms and a 
different one and another model where there are three different species of atoms. In both 
cases the three atoms can be combined to form a molecule. We then introduce a model 
consisting of three identical species of atoms which uses a recently defined multibosonic 
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realization of the si (2) algebra [IB]. We present how these generalized models can 
be derived from a transfer matrix, thus allowing us to show their integrability by the 
algebraic Bethe ansatz method. We also obtain their corresponding energies and Bethe 
ansatz equations. 

2. Triatomic molecular models 

We are going to introduce in this section three different Hamiltonians describing 
triatomic-molecular BECs, two for the hetero-atomic case and one for the homo-atomic 
case. We are considering the coupling parameters real, such that the Hamiltonians are 
hermitian. The U parameters describe the S'-wave scattering, the \x parameters are the 
externals potentials and Q is the amplitude for interconversion of atoms and molecules. 
The operators Nj, j = a, b, c, d are the number operators acting in the Fock space. 

2.1. Heteronuclear molecular models 

We can construct two models for hetero-atomic molecular BECs, one with two identical 
species of atoms and a different one and another model where there are three different 
species of atoms. 

2.1.1. Model for 2 atoms a and 1 atom b 

The Hamiltonian that describes the interconversion of a heterogeneous tri-atomic 
molecule labelled by c with two atoms of type a and one atom of type b is given by 

H = U aa N 2 a + U bb N 2 b + U CC N 2 C + U ab N a N b + U ac N a N c + U bc N b N c 

+ fi a N a + fi b N b + /j c N c + VL(a)a)tfc + c^baa). (2) 

This Hamiltonian (j2J) has two independent conserved quantities, 
h=N a + 2N C , 

h = N a - 2N b} (3) 

where I2 is the imbalance between the number of atoms of type a and b. The total 
number of particles N = N a + N b + 3N C can be written using these conserved quantities, 

2 ' 

and it is also conserved. 

Writing the S'-wave diagonal part of (j2D as a combination of the conserved quantities 
© we find 



alt + pil + jhh, 
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where we have used the following identification for the coupling constants 
U aa = a + (3 + 7, Uu = 4/3, U cc = 4a, 

U ab = -4(3 - 2j, U ac = 4a + 2 7 , U bc = -4 T . 
Therefore the Hamiltonian (j2j) can also be written as 

H = all + (3ll + 7/1 1 2 

+ fi a N a + /j b N b + /j c N c + fi(a t a t 6 t c + c ] baa). (4) 

This form will be used later when we show how the Hamiltonian ([2]) can be derived from 
a transfer matrix to establish integrability in the general context of the Yang-Baxter 
algebra. 

2.1.2. Model for 3 different atoms: a, b and c 

The Hamiltonian for a model of heterogeneous tri-atomic molecule labelled by d with 
one atom of type a, one atom of type b and one atom of type c is given by 



H = U aa N 2 a + U bb Nl + U CC N 2 C + U dd N 2 d 

+ U ab N a N b + U ac N a N c + U ad N a N d + U bc N b N c + U bd N b N d + U cd N c N d 
+ fi a N a + fi b N b + /j c N c + /j d N d + n(aW(2 + d ] cba). (5) 

This Hamiltonian (J5j) has three independent conserved quantities, 



h=N a + N d , 

I 2 = N b + N d , (6) 
h = N c + N d . 

The physical quantities representing the different imbalances between the number of 
atoms of different species (J ab , J ac , Jbc) and the total number of atoms N can be expressed 
as a combination of these quantities, for example, 



J ab = N a -N b = h- I 2 , 
J ac = N a -N c = h- I 3 , 
J bc =N b -N c = I 2 - I 3 , 
N =N a + N b + N c + 3N d = h + I 2 + I 3} 

and are conserved as well. 

We can write the S'-wave diagonal part of the Hamiltonian (JSj) in terms of a 
combination of the independent conserved quantities ([6]) as 
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all + (3I 2 2 + 61! + 7/1/2 + phh + OI 2 h, (7) 
where the following identification has been made for the coupling constants 

U aa = a, U bb = (3, U cc = 5, U dd = a + f3 + 5 + -f + p + 9, 

U ab = ~?, U ac = p, U bc = 9, tf M = 2/3 + 7 + 0, U cd = 26 + p + 9. 

Then the Hamiltonian (jSJ) can be also written as 

H = all + P 1 ! + 5I ! + 7/1/2 + phh + Ohh 

+ p a N a + p b N b + /i c A?" c + p d N d + J2(o t 6 + c t d + Scba). (8) 

Homonuclear-molecular model 

In the homogeneous case, we can construct a model describing a triatomic-molecular 
BEC, where all atoms are identical. 

2.2.1. Model for 3 atoms a 

The Hamiltonian that describes the interconversion of a homogeneous triatomic molecule 
labelled by b with three atoms of type a is given by 



H = U aa N 2 a + U bb Nl + U ab N a N b + p a N a + p b N b 

+ Q(a j a j a j ba-{N a ) + a„(N a )tfaaa). (9) 

where a-(N a ) is a function of N a (see next section for more details) that controls the 
amplitude of interconversion Q. This indicates that the density of atoms N a has some 
influence in the generation of a bound-state composed by three identical atoms. 
This Hamiltonian (Q has one conserved quantity, 



I = N a + 3N b , (10) 
the total number of particles N — N a + 3N b . 

3. Integrability and exact Bethe ansatz solution 

In this section we will discuss the derivation, the integrability and the Bethe ansatz 
solution of these models. We begin with the su(2)-invariant i?-matrix, depending on 
the spectral parameter u: 

( 1 \ 

b(u) c(u) 



R(u) 



c{u) b{u) 

y oiy 
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with b(u) = u/(u + rj) and c(u) = r\ j '{u + 77) . Above, 77 is an arbitrary parameter, to be 
chosen later. It is easy to check that R(u) satisfies the Yang-Baxter equation 

#12 (it - v)R 13 (u)R 23 (v) = R 23 {v)R n {u)R 12 {u - v). 

Here Rjk(u) denotes the matrix acting non-trivially on the j-th and the k-th spaces and 
as the identity on the remaining space. 

Next we define the Yang-Baxter algebra T(u), 



T(u) 



A(u) B(u) 
C(u) D(u) 



subject to the constraint 

R 12 {u - «)Ti(u)T 2 (t;) = Tz^T^Ruiu - v). (11) 

In what follows we will choose different realizations for the monodromy matrix 
7t(T(m)) = L(u) to obtain triatomic models for heteronuclear and homonuclear molecular 
BECs. In this construction, the Lax operators L(u) have to satisfy the relation 

R 12 {u - v)L x {u)L 2 {v) = L 2 {v)L 1 {u)R 12 {u - v). (12) 

Then, defining the transfer matrix, as usual, through 

t( u ) = tm(T(u)) = n{A(u) + D(u)), (13) 

it follows from ffTTl) that the transfer matrix commutes for different values of the spectral 
parameters; i. e., [t(u),t(v)] = 0. Consequently, the models derived from this transfer 
matrix will be integrable. Let us now particularize this construction for the hetero and 
homo-atomic molecular cases. 

3.1. Heteronuclear-molecular models 

We may choose the following realization for the two models of heteroatomic molecular 
BECs 

tt(T(u)) = L{u) = u-GL s (u~)L K {u + ), (14) 
with w ± = u ± u, G = diag(+, — ), the Lax operator L (u) 

u \ —r/b u + rjb 

in terms of the su(2) Lie algebra with generators S z and S^ 1 subject to the commutation 
relations 

[^ 2 ,S ± ] = ±S ± , [S + ,S'} = 2S Z , 
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and the Lax operator L K 
L K (u) 



K 1 / U + T]K Z 7]K 



u \ —r]K + u — r/K z 

in terms of the su(l, 1) Lie algebra with generators K z and K ± subject to the 
commutation relations 

[K Z ,K ± ] = ±K ± , [K + , K-] = -2K Z . (15) 

Now, using two different realizations for the su{2) and su(l, 1) algebras we will show 
how to construct the hetero-atomic Hamiltonians ([2]) and ([5]) from the transfer matrix 
( JT3l) and present their exact Bethe ansatz solution. 



3.1.1. Model for 2 atoms a and 1 atom b 

Using the following realizations for the six (2) and su(l, 1) algebras 

N h -N r . 



S+ = 6+c, S~ = c\ S z 



2 



2 ' 2 ' 4 ' 

it is straightforward to check that the Hamiltonian is related with the transfer matrix 
t{u) (P3J through 

H = t(u) - ^u-f] + all + (3ll + 5/i/ 2 , 
where the following identification has been made for the parameters 

fi a = u~r], fi c = -fi b = u + 7], Q = — . 

We can apply the algebraic Bethe ansatz method, using as the pseudovacuum the 
product state (|0) = \k) <S> \<j>), with \k) denoting the lowest weight state of the su(l, 1) 
algebra with weight k, i.e., K z \k) = k\k) and \(f>) denoting the highest weight state of 
the sit (2) algebra with weight m z ), to find the Bethe ansatz equations (BAE) 

Jp-,- V n,) { v, + . + V k) ^v 1 ~v^^ { v 

[Vi — uo + r]m z ) [Vi + u> — rjk) -V A Vi — Vj + r\ 

and the energies of the Hamiltonian (j2J) (see for example [7J [8] ) 

M u _ v 

E = (m _ — r]m z ){u + + i]k) J^J 

i=i 

M 

+ r]m z )(u + — rjk) 



U-Vi 

i=i 



M 

u — Vi — ri 



U-Vi 

i=i 



- -vT-n + all + pi! + 7/1/2. (17) 
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3.1.2. Model for 3 different atoms: a, b and c 

Using the following realizations for the su(2) and su(l, 1) algebras 

N — 

S+ = c% S- = dfic, S z = c d , 

K+ _ (aW) ir= W K ,_ N a + N b + i 

2 ' 2 ' 2 

it is straightforward to check that the Hamiltonian (jSJ) is related with the transfer matrix 
(USD through 

H(u) = t(u) - iTrj + all + f^l + <"! + 7/1/2 + p/1/3 + 0/2/3, 
with the following identification for the parameters 

Applying the algebraic Bethe ansatz method we find the Bethe ansatz equations 
(BAE) 

(Vi - u - rim z )(vi + u + -qk) ^Vj-Vj-y . , . 

-7 77-— m = 11 ; — > «>J = (18) 

and the energies of the Hamiltonian (JSJ) 

M M - U- + 

i£ = (it - — rjm z )(u + + 77/c) I I 

"7 w - f » 

M 

u — Vi — ri 



+ rjm z )(u + - ink) JJ 



- vTr) + all + (3I 2 2 + 5/f + 7/1/2 + Phh + Ohh- (19) 
3.2. Homonuclear-molecular model 

In this case we may choose the following realization for the Yang-Baxter algebra 

tt(T(u)) = L{u) = rj- 1 GL b (u -5- r]- l )L A (u + u), (20) 
with, G = diag(—, +), the Lax operator L b (u) 



L b (u) 



u + rjN b b 



6t 



in terms of the canonical boson operators b and with N b = b'b, subject to the 
commutation relations 



[b,b] = [tf,tf]=0, [&,&+] = !, 
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and the Lax operator L A 

L A (u)=( U + l A ° 

\ -r)A + u-^Ao ' 

in terms of the si (2) Lie algebra with generators A and A±, subject to the commutation 
relations 

[A., A + ) = Ao, [A , A±] = ±2A±. (21) 
We find from the L operator (I20I) the following transfer matrix 

t( u ) = - rrV + w + ~A )(u - 5 - r]- 1 + 7]N b ) 
V 

■f r/ - \u -t- — 

with 



+ r7- 2 (M + cu-^ ) + M + + 6U_, (22) 



t(0) = ^-^(5 + 2r]^) - cuiV fe + ~(5 - r)N b )A + bA + + (23) 



There is a one-mode realization of the si (2) algebra 

A = a (N), A_=a4N)a l , A + = (a^) l a„(N) 



with 

ao(JV) =y(iV- J R) + « ( J R) (24) 
/ iVI 1 1 

a - (iV) = V (ivT7)I ( T (iV ~ H) + ao(jR))( T (iV " ^ + 1} (25) 

where N = a^a and / G N. The operator R is 

0, if 1 = 1 

~2 ^ e (27rm/0 _ X ' ^ 

m=l 

and acts on the states {|n)} as R\n) = n mod l\n). The function ao(R) is a function of 
the spectrum of R. For n = r < I, we have 

j{N-R)\r) = 0|r) (26) 

with A = «o(-R) such that ao(R)\r) = ao(r)|r) and R\r) = r\r). The operator R 
commutes with a 1 and (a^) 1 and so is a conserved quantity for all models presented in 
this paper. 

Now we will apply this realization to show how to construct the Hamiltonian 
from the transfer matrix (1221) and present their exact Bethe ansatz solution. For I = 3, 
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which is the case of interest here, it is straightforward to check that the Hamiltonian 
(Q is related with the transfer matrix t(0) fl23l) through 

H = t(0) (27) 

where we have the following identification 



v = 9U a + U b - 3U t 



ab 



6(u + 5) = (18U a - 2U b )N -r) P + 18/i a - 6// 6 

36^ _1 a;5 + 72w»T 2 = r]p 2 + [{AU b - 6U ab )N + 6(p b - 3p a + uj)] P 
+ 4U b N 2 + 12^ + oo)N 

with p = p(-R) = 3ao(-R) — 2_R. Let |0)& denotes the Fock vacuum state and let 
\t)a denotes the lowest weight state of the si (2) algebra where r = 0, 1, 2, are the 
eingenvalues of R for I = 3 and N = nl + r, with n G iV. On the product state 
1^) — l r )yi <8> |0)b we can a Pply the algebraic Bethe ansatz method to find the Bethe 
ansatz equations (BAE) 

(1 - rjVi + r)6)(vi + u + %a (r)) ^L Vi - Vj -r] . . 

; * — r~\ = ll r~> hJ = h-,M, (28) 

and the energies of the Hamiltonian (j2j) 



M 



E = rf\8 + rT 1 ) (w + |a (r)) JJ 



+ .r 2 (--%(r))n^- (29) 

1=1 

4. Summary 

We have introduced three new integrable models for both, homogeneous and 
heterogeneous tri-atomic molecular BECs obtained through a combination of Lax 
operators constructed using special realizations of the su{2) and su(l, 1) algebras 
and a particular one-mode multibosonic representation of si (2), possibilities that were 
overlooked in previous studies. The models were solved by means of the algebraic 
Bethe ansatz method and their corresponding energies and Bethe ansatz equations were 
derived. 
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